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A general covariant extension of Einstein's field equations 
is considered with a view to Numerical Relativity applica- 
tions. The basic variables are taken to be the metric tensor 
and an additional four-vector. The extended field equations, 
when supplemented by suitable coordinate conditions, deter- 
mine the time evolution of all these variables without any con- 
straint. Einstein's solutions are recovered when the additional 
four-vector vanishes, so that the energy and momentum con- 
straints hold true. The extended system is well posed when 
using the natural extension of either harmonic coordinates or 
the harmonic slicing condition in normal coordinates 



I. INTRODUCTION 

Einstein's field equations can be understood as a set of 
ten second order partial differential equations on the ten 
unknown metric coefficients g^i, 
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Rtiv = 8 TT (T^i, - -T g^^) 



(1) 



But in General Relativity, like in other field theories, 
physical solutions are given instead by an equivalence 
class of mathematical expressions which are related one 
to another by gauge transformations. The peculiarity of 
General Relativity is that the gauge group is that of the 
general (smooth) coordinate transformations 



(2) 



This allows one to use the coordinate gauge freedom (g) 
to fix up to four of the ten metric coefficients (kinematical 
degrees of freedom) so that only the remaining ones (dy- 
namical degrees of freedom) can be actually determined 
by the field equations (|l|). It follows that four of these 
equations just provide additional constraints on the dy- 
namical degrees of freedom: they can be understood as 
first integrals related with energy and momentum con- 
servation. 

This complex structure becomes more transparent in 
the context of the 3-1-1 formalism where one consid- 
ers the space-time sliced by t = constant hypersurfaces. 
The line element can then be written as 



-a^ dt^ 



+ "f^j {dx' + dt) {dx^ + 13^ dt) i, j = 1, 2, 3 (3) 



where the lapse a and the shift /3* are the kinematical 
degrees of freedom, whereas the coefficients "fij of the 



metric induced on the t ~ constant hypersurfaces are 
the dynamical ones. One can use the extrinsic curvature 
(second fundamental form) of these slices 
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Kij ^ -—{dt - Cp) -fij 



(4) 



to decompose the field equations (^ into six evolution 
equations 



r(3) 



2Kt 



(ADM evolution equations) and four constraints 

^^'^R-tr{K^) + {trKf =0 
Vfe (if^ - trK ,5^) = 



(5) 



(6) 
(7) 



where we have restricted ourselves to the vacuum case 
for simplicity. 

The 3-1-1 formalism (^-0) is specially suited for Numer- 
ical Relativity applications. In this context, one usually 
takes advantage of the fact that the energy and momen- 
tum constraints (^J^) are first integrals of (^,||). This 
allows one to impose (^,^ on the initial data only ('free 
evolution' approach), but not during time evolution. One 
is dealing then with an extended set of solutions. The 
resulting simulations will represent true Einstein's solu- 
tions only to the extent to which the constraints (HJ^) are 
actually preserved during numerical evolution. 

This approach raises some concerns (seeRef. [|). First 
of all, the set of 'extended solutions' is not determined 
in a general covariant way: equations (^,^ just translate 
the space components of the original field equations (j^). 
And the space components of four-dimensional tensors 
are covariant only under the restricted subgroup of co- 
ordinate transformations that preserve the time slicing 
(3-1-1 covariance): 



t'^h{t), f = f{x^,t) 



(8) 



Moreover, even in a given coordinate system, there are 
many ways of extending the space of solutions, because 
the set of six evolution equations (||) is not univocally 
defined in that context: one can actually combine either 
the energy constraint (^ or the momentum constraint (|^) 
with the original evolution equations (|^) to get a different 
set of evolution equations. This fact has been recently 
used by Kidder, Scheel and Teukolsky (KST) to obtain 
a multiparameter family of evolution systems [Q. 
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This ambiguity was previously used in a different way 
by Bona and Mass [^,D, and also by Shibata and Naka- 
mura and Baumgarte and Shapiro (BSSN system). 
The key idea in these works was to introduce three sup- 
plementary dynamical variables whose evolution equa- 
tions were obtained by using the momentum constraint 
. As far as these works were focused on Numerical Rel- 
ativity applications, the supplementary quantities were 
introduced in an 'ad hoc' way, breaking the 3-1-1 covari- 
ance of the formalism. Only very recently |^ the same 
idea has been implemented in a 3-f 1 covariant way: the 
extra quantities are given by a three-dimensional 'zero' 
vector Zi which vanishes for Einstein's solutions. During 
numerical evolution, however, non-zero values of Zi arise 
due to truncation errors so that the resulting numerical 
codes actually deal with an extended set of solutions. 

The next logical step along this line would be to in- 
troduce one more supplementary quantity, let us call it 
0, whose evolution equation would be obtained by using 
the energy constraint in the same way as momen- 
tum constraint was used in |9| to evolve Zi. The new 
quantity could be then coupled to the other equations 
in many different ways so that many more new parame- 
ters will appear in the resulting evolution system (see for 
instance Ref. jl^] for a similar approach). 

We will not go that way in the present work. We prefer 
to preserve general covariance in our extension of Ein- 
stein's field equations. Instead of a separate three-vector 
Zi and three-scalar 0, we shall consider a covariant 'zero' 
four- vector such that their space components coincide 
with Zi and 



where we have rescaled the four-vector to clean ( [Tl| ) 
from any arbitrary parameter. The Cauchy problem of 
the extended equations (^ij) is simpler than the corre- 
sponding one for (|l|). The space components of ( pT| ) 
provide again second order evolution equations for the 
dynamical degrees of freedom in the metric. In the 3-1-1 
language, one has 

idt - Cp)K,j = -V,aj + a f^R^j + V,Z, + V^Z, (12) 
- 2 Kl + {trK - 2 6) K,^] 

But now the remaining components can be combined 
with (^2|) to provide, instead of constraints, first order 
evolution equations for the components of Z^, namely 

{dt- Cp)Q='^[^^'iR+{trK -2Q)trK (13) 

- ir(i^2) + 2 VkZ^ - 2 Z^akla] 

{dt - Cp)Z, = a [V, (K,^ - S.HrK) + 9,9 (14) 

- 2 Ki^ Zj - 6 aja] 

where here again we have restricted ourselves to the vac- 
uum case. 

The contracted Bianchi identities, when applied to 
eq.(|l^), lead to a simple equation for 



(15) 



where stands for the four-dimensional Lichnerowicz 
Laplacian. This equation, of course, can be also derived 



e = n„ Z'' = a Z° 



(9) 



from (|12-14), but the converse is not true: the first order 
evolution equations (^ 14) for Z^ can not be derived 
from ([l^) and the second order equation (p^). 



where is the unit normal to the t = constant slices. 
The purpose of this work is to use Z^ to extend in a 
covariant way Einstein's equations so that the extended 
system contains only evolution equations, without any 
constraint. Einstein's solutions will be recovered for 
Z^ = 0. The four- vector Z^ will then provide a simple 
way to measure of the quality of numerical simulations: 
energy constraint violations can be monitored with O^, 
whereas one can use the norm of Z^ to deal with mo- 
mentum constraint deviations. These quantities can be 
combined to form a four-dimcnsional scalar: 



z^z'' = - 7*-'z,Zj 



(10) 



III. COORDINATE CONDITIONS 



In order to analyze the causal structure of the extended 



equations (11 ), we will use de DcDondcr ||ll|] expression of 
the Ricci tensor to write down the principal part, namely 



□ 



Sm-^ + ^^(r. + 2 z,) + a,(r^ + 2 z^) 



(16) 



where here the box symbol stands for the d'Alembert 
operator on functions and we have noted as usual = 
gP'^Y^p^. Comparing with the wave equation for 17^1^, it 
is clear that we can obtain a well posed system if we take 
out the additional terms in ( [l6| ) by using the following 
extension of the well known harmonic coordinate condi- 
tions: 



II. EXTENDING THE FIELD EQUATIONS 



-r^ = 2 z'" 



(17) 



We propose to extend the field equations (|l|) in the 
following way 

i?,,. + V,,Z, + V,Z^ = 8 TT (T^, - ir g^,) (11) 



so that the four-vector Z^ can be interpreted in this con- 
text as providing 'gauge sources', along the lines sketched 
in ref |l|. 

Harmonic coordinates, including their extension (p^), 
are not flexible enough to be used in most Numerical 
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Relativity applications. A more suitable choice is the 
'harmonic slicing', in which the time coordinate is again 
assumed to be an harmonic function, but the space co- 
ordinates are chosen so that the mixed components goi 
vanish. We propose here to keep in the extended case the 
time component of (p7|), that is 



ax" 



2Z" 







(18) 



which can be translated into the 3+1 language as 
{dt - Cp) \na = ~a [trK - 2 9), = 



(19) 



The mixed time-space components of (16) read now 



(20) 



This means that in the harmonic slicing case (18) there 
are three degrees of freedom, directly related with the 
space components of Z , with characteristic lines orthog- 
onal to the t = constant slices. This 'standing' modes 
can be decoupled from the remaining ones by taking an 
extra time derivative of (|l^). This means that the re- 
sulting system will be of third order in a and gij (second 
order in 8 and Kij). Allowing for ( p^ , p^ , one gets after 
some algebra 



□ 6 = 



(21) 



so that the corresponding characteristic lines are on the 
light cones. It follows from (|2^j2^) that the resulting 
system (with that extra time derivative) is well posed (see 
Ref. |l^ for a similar development, in the 3-1-1 language, 
from Einstein's equations (Q)). 

Although numerical simulations are beyond the scope 
of the present letter, we will comment here some exten- 
sions of this work which are specially suited for Numer- 
ical Relativity. The first one is that harmonic slicing is 
often generalized in that context by including an extra 
factor /(a) in the right-hand-side of equation (p^. The 
effect of this minor change is twofold: from the theoreti- 
cal point of view, the gauge-related degrees of freedom no 
longer propagate with light speed (gauge speed) so that 
the first equation in ( pl| ) no longer holds true; from the 
numerical point of view, higher gauge speed increases the 
singularity avoidance properties of the gauge |l^] . 

The second extension is to consider the first space 
derivatives 



to Numerical Relativity applications than the above men- 
tioned third order system. In addition, the causal struc- 
ture of first order systems is much easier to analyze |^ 
and one can take advantage of the characteristic field de- 
composition to devise stable numerical algorithms. 
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as independent dynamical variables with evolution equa- 
tions 



dtAk + dk[a f{trK ~2Q)]=Q 
dtDkij + dk[a. K,j] = 



(23) 
(24) 



so that the set mM 



24 ) forms a first order hyperbolic 



evolution system, which is obviously much more adapted 
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